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This proceeding reports a study of the impacts of the SB2009 beam parameters on the 
Higgs recoil mass and Higgs-Strahlung cross-section measurements based on a dedicated 
fl_) ■ fast simulation algorithm of the ILD detector. The study shows worse results from 

SB2009 beam parameters than the previous RDR beam parameters, because of the 
smaller luminosity. However, Travel Focus (TF) technology can recover the degradation 
to certain level. 



1 Introduction 

>, 

The intent of this work is to study of the impacts of the SB2009 [TT] beam parameters on 
the Higgs recoil mass and Higgs-Strahlung cross-section measurements based on a dedicated 
fast simulation algorithm of the ILD [2] detector. We would like to compare this result with 
our previous result [1 H S] at y/s =250 GeV with RDR [1] beam parameters (RDR 250) for 
the accelerator optimization reason. 

The study uses the ZH — > (i + fi~X channel, assuming a beam polarization of (e~ : 
—80%, e+ : +30%). It takes into account the beam effects by beam simulation using 
GUINEA-PIG [8] with beam parameters SB2009. I use PYTHIA [S] for the event generation 
with the beams simulated by GUINEA-PIG as inputs through the interface CALYPSO [TU] . 
Thereafter, I developed a fast simulation algorithm [5] to include the detector effects of the 
ILD. After the fast simulation, I perform the analysis based on the same algorithm as of my 
previous study J21EJ- In the end, I give the result and the comparison. 

2 Beam Simulation 

The beam effects are simulated using GUINEA-PIG, with the beam parameters SB2009 
given by reference |llj . Figure [1] shows the luminosity spectra simulated with beam param- 
eters SB2009, comparing with RDR beam parameters. . 

3 Event Generation 

I use PYTHIA for the event generation. The event generation takes the beams simulated 
by GUINEA-PIG as inputs through the interface CALYPSO. 
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Figure 1: Normalized luminosity spectra simulated with beam parameters of SB2009 and 
RDR. 

I take the ZH — > /i + [i~X channel under study, with two major background reactions 
the WW (W + W~ ->■ fi+v^/jTiTjj,) and the ZZ {ZZ -+ (J, + /J,~ff). Their cross-sections at 
t/s =350 GeV with beam polarization (e~ : —80%, e + : +30%) are listed in Table [TJ 



Reaction 


Cross-Section 


ZH ->• /jftX 


7.1 fb 


WW 
ZZ 


346 fb 
165 fb 



Table 1: Reactions and cross sections at y/s =350 GeV with beam polarization (e : 
—80%, e + : +30%). The signal is indicated by bold face letters. 

If I take the RDR 500 peak luminosity (£ pea k,RDR500 = 2. Ox 10 34 cm _2 s _1 ) and integrated 
luminosity (Ant,RDR500 =500 fb -1 ) as reference, the estimated integrated luminosity of a 
given set of beam parameters should be [TJ] : 



£int — 



£ 



peak 



c 



■ c 



poak,RDR500 



int,RDR500 



(1) 



Following this rule, these integrated luminosities for various beam parameters are listed 
in Table [U 

4 Fast Simulation 

I developed a dedicated fast simulation algorithm for the ILD detector concept. 
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RDR 


SB2009 w/o TF 


SB2009 w/ TF 


v^ (GeV) 


250 


350 


500 


250.a 


250.b 


350 


500 


250.a 


250.b 


350 


500 


Peak L (lO^cm"^- 1 ) 


0.75 


1.2 


2.0 


0.2 


0.22 


0.7 


1.5 


0.25 


0.27 


1.0 


2.0 


Integrated L (fb _I ) 


188 


300 


500 


50 


55 


175 


375 


63 


68 


250 


500 



Tabic 2: Estimated Integrated luminosities for various beam parameters |11| . 



The fast simulation of the ILD detector intents to consider the uncertainty of the detector 
response without full detector simulation. For the Higgs recoil mass analysis, the major 
detector uncertainty is due to the momentum measurement of the lepton tracks. I take the 
muons tracks for this study. Therefore the object of the fast simulation is to smear the MC 
true momentum given by event generation according to the momentum resolution of the 
detector. It consists of two steps. 

The first step is to have the momentum resolution of the ILD detector. I parameterize 
the momentum resolution of the ILD detector as a function of the momentum (P) and cos 9 
of muon tracks. The momentum resolution function is given by Equation^ 



I oi©6i/P : |cos0| < 0.78 

P2" 1 (a 2 e&2AP)/sm(l-|cos6>|) : |cos6»| > 0.78 



(2) 
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Figure 2: The black scatter plot shows the distribution of momentum resolution of the ILD 
detector as a function of lepton P and cos 9 obtained from full simulation of the detector. 
The red surface shows the fit of Equation [2] to the distribution. The parameters obtained 
from the fit is shown in Table El 

Equation [5] is thus fitted to the distribution of momentum resolution got from the full 
simulation of the ILD detector, as shown in Figure [5] In this figure, the black scatter 
plot shows the momentum resolution as a function of P and cos 9 from ILD detector full 
simulation, and the red surface is the fit of Equation [2] to the distribution of momentum 
resolution. The momentum resolution for a given lepton is obtained by propagation of the 
error matrix of its track reconstruction. The fitted parameters are shown in Table [3] 
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dl 


2.08 x l(T a (1/GcV) 


h 


8.86 x 10- 4 


02 


3.16 x 10~ b (1/GeV) 


02 


2.45 x 1(T 4 



Table 3: Parameters in Equation [2] obtained by fitting it to the distribution of momentum 
resolution of the ILD detector shown in Figure [2j 



The second step of the fast simulation is to smear the MC true momentum from the 
event generation according to the momentum resolution function just obtained. For each 
lepton, I define a Gaussian function, and set the MC true momentum as the mean of this 
Gaussian function, the AP as the sigma. Thereafter, I generate a random number according 
to this Gaussian function, and this random number is the momentum after fast simulation. 
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Figure 3: Comparison of the invariant mass Mdi of the lepton pair (top) and the recoil mass 
Mh (bottom) distributions from fast simulation and full simulation of the ILD detector at 
t/s =250 GeV. 

In order to validate the method, I compare the Mdi and Mh distributions of this fast 
simulation with the full simulation at -y^ =250 GeV. Figure [3] shows this comparison. They 
agree with each other. I have also done a further validation by repeating the analysis at 
■s/s =250 GeV for the ILD LOI based on this fast simulation, and comparing the results 
with those from the full simulation. The results are identical to each other. 

Figure 0] shows the comparison of the Higgs recoil mass distributions of the signal af- 
ter fast simulation with various beam parameters considering the corresponding integrated 
luminosity. 

5 Analysis and Results 

The analysis procedure [3] after the fast simulation is similar to that of my study based on 
the full simulation samples: 

• A cut based background suppression. The cuts are defined in Table 2] These cuts are 
independent of the Higgs decay mode. 
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Figure 4: Comparison of the Higgs recoil mass distributions of the signal after fast simulation 
with various beam parameters considering the corresponding integrated luminosity. The 



beam polarization is (e 



-80%, 



-30%: 



• A Likelihood further rejection of background j3]. The variables (Prdl, cos 6 di, Mdi and 
acolinearity ) employed in this Likelihood suppression are also independent of the 
Higgs decay mode. 

• The resulting recoil mass (M reco u) spectrum of signal and background is fitted to 
derive the results of the Mh and the ZH cross-section measurement [3] . I choose the 
physics motivated function [3] to describe the signal. 



Cut-Chain 



(1) |cos0 M | <0.99 

(2) P Tdl > 20 GeV 

(3) M dl e (80, 100) GeV 

(4) acop e (0.2, 3.0) 

(8) M recoil G (115, 150) GeV 

(9) Likelihood Further Rejection 

(using variables Prai-, cos 9 di, Mdi and acol) 

Table 4: The cut-chain for background suppression. 

After the background suppression, the remaining numbers of events of signal and back- 
ground reactions are given in Table [5] 

Figure [5] shows the fit to the recoil mass spectrum of remaining signal and background. 
An equivalent plot with beam parameters RDR 250 can be found in reference |31 0] . 

From Figure [5] the derived results of the Higgs mass measurement is Mh = 120.010 ± 
0.110 GeV, and of the cross-section is a — 7.13 ± 0.37 fb (5a ja = 5.2%). A comparison of 
the results with other beam parameters are shown in Table [6l together with the efficiency 
(e) and signal over background (S/B). In this table, the results of RDR 250 SB2009 w/o 
TF 250b are scaled from my previous analysis [3l H] based on full simulation according to 
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Reactions ZH -> \x\iX ZZ WW 



^initial 

■t-V selected 



1248 
633 



29k 
658 



61k 
30 



Table 5: Numbers of events before and after background suppression, for signal and back- 
grounds. 
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Figure 5: Fit to recoil mass spectrum of signal plus background based on the fast simulation, 
at y/s =350 GeV, with beam polarization (e~ : —80%, e + : +30%) beam parameters SB2009 
w/o TF, and assuming an integrated luminosity of 175 fb -1 . 



Beam Par 


Ant (fb- 1 ) 


e 


S/B 


M H (GeV) 


a (fb) (5a/ a) 


RDR 250 


188 


55% 


62% 


120.001 + 0.043 


11.63 + 0.45 (3.9%) 


RDR 350 


300 


51% 


92% 


120.010 + 0.084 


7.13 + 0.28 (4.0%) 


SB2009 w/o TF 250b 


55 


55% 


62% 


120.001 + 0.079 


11.63 + 0.83 (7.2%) 


SB2009 w/o TF 350 


175 


51% 


92% 


120.010+0.110 


7.13 + 0.37 (5.2%) 


SB2009 TF 250b 


68 


55% 


62% 


120.001 + 0.071 


11.63 + 0.75 (6.4%) 


SB2009 TF 350 


250 


51% 


92% 


120.010+0.092 


7.13 + 0.31 (4.3%) 



Table 6: Results of different beam parameters, assuming a beam polarization of (e : 
-80%, e+ : +30%). The results of RDR 250 and SB2009 w/o TF 250b are scaled from my 
previous analysis based on full simulation according to the integrated luminosity. That of 
RDR 350 is estimated by scaling the result of SB2009 w/o TF 350 obtained here according 
to the integrated luminosity. 
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the integrated luminosity. And the results of RDR 350 is estimated by scaling the result of 
SB2009 w/o TF 350 according to the integrated luminosity. 

The higher S/B at y/s =350 GeV is due to better background suppression. For example 
the variable cos 8di, its distribution of ZH signal is much center for y/s =350 GeV than 
250 GeV [3], while that of the ZZ background is much forward. 

When comparing RDR 250 and 350, the errors on the cross-section are similar, while the 
error on the Mh at 350 GeV is worse by a factor of 2 than that at 250 GeV. 

When comparing SB2009 w/o TF 250 and 350, the higher luminosity at 350 GeV gives 
better result on cross-section, but not on the Mh, which is worse by 1.4 times due to the 
wider mass peak. 

For a given s/s, the results of SB2009 w/o TF are worse due to the decrease of luminosity. 
The TF recovers the degradation to certain level. 

Also, the comparison above shows the results on Higgs mass is about 3 times worse if 
we change to use beam parameters SB2009 350 w/o TF, and on the cross-section more than 
1.5 times worse. 
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